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ABSTRACT 

We present the results of an unbiased radio search for gravitational lensing events 
with image separations between 15" and 60", which would be associated with clusters 
of galaxies with masses > 10^^ ~ ^^Af©- A parent population of 1023 extended radio 
sources stronger than 35 mJy with stellar optical identifications was selected using 
the FIRST radio catalogue at 1.4 GHz and the APM optical catalogue. The FIRST 
catalogue was then searched for companions to the parent sources stronger than 7 mJy 
and with separation in the range 15"to 60". Higher resolution observations of the 
resulting 38 lens candidates were made with the VLA at 1.4 GHz and 5 GHz, and 
with MERLIN at 5 GHz in order to test the lens hypothesis in each case. None of our 
targets was found to be a gravitational lens system. These results provide the best 
current constraint on the lensing rate for this angular scale, but improved calculations 
of lensing rates from realistic simulations of the clustering of matter on the relevant 
scales are required before cosmologically significant constraints can be derived from 
this null result. We now have an efficient, tested observational strategy with which 
it will be possible to make an order-of-magnitude larger unbiased search in the near 
future. 

Key words: cosmology: observations, (cosmology:) gravitational lensing, galaxies: 
clusters: general, radio continuum: general 



1 INTRODUCTION 

Observations of gravitational lensing provide a probe of the 
distribution of mass in the universe on scales ranging from 
micro-lensing by sub-solar mass objects in the Galactic halo 
to weak-lensing by large-scale structures of the universe. In 
order to survey matter both light and dark on different size 
scales, unbiased systematic searches for lensing events within 
a well-defined sample of background sources must be carried 
out. Since the discovery of giant arcs in the region of clus- 
ters A370 and C12244 (Lynds & Petrosian 1986; Soucail et 
al. 1987) more than 50 clusters have been found to contain 
arc-like images of background galaxies, including a clear ex- 
ample of multiple imaging by the cluster 0024-f 1654 (Col- 
ley, Tyson & Turner 1996). The lensed images provide di- 
agnostics for the cluster potential, allowing the distribution 
of mass to be mapped independently of the distribution of 
light. However, the currently known examples cannot tell us 
about the general prevalence of mass structures on the scale 
of clusters of galaxies (CGs) because they are in no sense sta- 
tistically complete, having largely been discovered through 
observations of previously known clusters. Only by carefully 
examining a large well-defined sample of background sources 
for cases of gravitational lensing can we hope to place con- 
straints on the general distribution of cluster-sized masses. 



A''-body cold-dark-matter (CDM) simulations of the 
universe make predictions about the degree of clustering of 
dark matter haloes at different epochs (for example Jenk- 
ins et al. 1998). The statistics of strong gravitational lens- 
ing events within a sample of background sources should 
therefore allow us to distinguish between models with dif- 
ferent rates of cluster evolution. Early simulations by Wamb- 
sganss et al (1995) suggested that there could be as many 
strongly-lensed sources with image separations >5", associ- 
ated with groups and clusters of galaxies, as there are with 
image separations <5" associated with single galaxy lenses. 
The JVAS and CLASS lens surveys (King et al. 1999; My- 
ers 1996) found a lensing rate of approximately 1 in 500 for 
lens systems with image separations between O'.'S and 6". 
This suggests that in a search for lens systems associated 
with clusters of galaxies, a background source population of 
at least ~ 10'^ is necessary before useful cosmological con- 
straints can be set. 

Previous searches for strong lensing events with image 
separations greater than a few arcseconds have started with 
quasar surveys because they provide a sample of objects 
with high mean redshift, which increases the lensing proba- 
bility, and relatively low surface density, which reduces the 
number of false positives from random associations. Maoz 
et al. (1997) used a subsample of 76 quasars from the HST 
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snapshot survey (Maoz et al. 1992) to look for lensing events 
with angular separations of 7" to 50". The fact that no lens 
systems were found allowed some extreme cosmological sce- 
narios to be ruled out, for example that there exists a large 
population of 'failed' or 'dark' clusters. However, the small 
size of the parent sample limited the ability to place useful 
constraints on more conventional cosmological models. 

In designing surveys to look for cluster-mass lenses, 
a basic problem arises when trying to recognise strongly- 
lensed images. As the image separation increases, the time 
delay between the images increases, and at a separation of 
tens of arcsec the characteristic delay can be hundreds of 
years. Basic properties of quasars such as the continuum 
emission, the emission line spectrum and the radio source 
nuclear structure may change significantly over this time. 
Due to their intrinsic variability, quasars are therefore not 
optimally-suited for use as probes for lensing by clusters of 
galaxies. A search method using background sources which 
are expected, a priori, to show no significant changes in phys- 
ical properties over timescales of ^1000 years is required to 
overcome this problem. 

This paper is concerned with the Arcminute Radio 
Cluster-lens Search (ARCS), a project aimed at looking for 
strong gravitational lensing by galaxy clusters in a system- 
atic and uniform manner. We use the fact that extended 
emission (>few kpc) from an extragalactic radio source is 
unlikely to vary over the timescales of interest. In order to 
maximise our parent sample of background sources, whilst 
minimising the observational effort in following up potential 
lens candidates, we started from the existing FIRST 1.4 GHz 
survey (Becker, White & Helfand 1995) and the Automatic 
Plate Measuring (APM) optical catalogue (for example Ir- 
win, Maddox & McMahon 1994). The FIRST survey pro- 
vides both high sensitivity and sufficient resolution to allow 
us to reject most multiple sources as candidate lens sys- 
tems directly before any follow-up observations have taken 
place. Furthermore, demanding stellar optical identifications 
of the radio sources from the APM catalogue both increases 
the mean redshift of the parent sample and provides a filter 
against selecting intrinsic double sources as lens candidates. 

We describe our selection algorithm for identifying grav- 
itational lens candidates from the FIRST and APM cata- 
logues in section g. We then go on to describe the subsequent 
follow-up observations in section 0, including an analysis 
of the remaiinng CG-lens candidates at each observational 
stage. A review of the project and future work is laid out in 
section n. 



2 LENS CANDIDATE SELECTION 

Our aim is to maximise the probability of detecting grav- 
itational multiple images produced by CGs while minimis- 
ing the observational effort required to confirm or refute 
them. As we pointed out above, reliable searches for lensing 
events on these scales require us to select background parent 
sources which should not vary over the relevant timescales 
which could extend to hundreds of years. We must also 
meiximise the redshift of these background sources which in- 
creases the lensing probability. The FIRST survey combined 
with the APM catalogue derived from the Palomar Observa- 





FIRST radio source 



1.5" window 



APM optical identification 



Figure 1. Illustration of a one-component (left) and two- 
component (right) primary source. Note that the position of the 
optical identification for a two-component source is determined 
by the relative brightness of each component (see texti . A range 
of two-component sources are acceptable (see Figure bh ranging 
from those with optical identification close to a relatively bright 
component, to those with two components of similar brightness 
and optical identification equidistant from each component. 




Centroid position 




Figure 2. Illustration of two extreme examples of two-component 
primary sources. The left-hand source has a centroid position near 
to a compact but relatively strong component whilst the second 
component is relatively extended but with significantly lower in- 
tegrated flux density. This type will most likely be a 'core-jet' 
source. The source on the right comprises two components with 
the same strength, and the centroid position is approximately 
equidistant from each component. In this case the source is prob- 
ably a 'classical double'. 



tory Sky Survey (POSS) plates enables us to achieve these 
goals. 



2.1 The FIRST Survey 

The Faint Images of the Radio Sky at Twenty centimetres 
(FIRST) survey will eventually cover 10 000 square degrees 
of the north galactic cap. The survey area will coincide with 
the Sloan Digital Sky Survey (SDSS), upon completion of 
which ~50 per cent of all FIRST sources will be optically 
identified to the m{v) ~ 24 limit of SDSS. For our pur- 
poses, the north Galactic cap portion of the catalogue dated 
98feb04 was used. This contains 382 892 sources covering 
4150 square degrees of the sky in the region RA~ 7*^20™ to 
RA~ n^2Cr, and 5 ~22'=12'to 5 ~57°36'. The lower flux 
density limit is 1 mjy with a typical rms noise of 0.15 mjy 
and a beam-size of 5'.'4; source positions are tied to the JVAS 
calibrator reference frame (Patnaik et al. 1992; Browne et al. 
1998; Wilkinson et al. 1998) with a systematic uncertainty 
of <0'.'03, while the overall systematic positional errors are 
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<0'.'l throughout the survey area. Source parameters were 
derived by fitting elliptical Gaussian models. A more com- 
plete description of the FIRST c atalogue can be found at 
URL ittp://sundog. stsci .edu/. 

Approximately 2.5 per cent of 'sources' in the FIRST 
catalogue are flagged as possible side-lobes of strong sources. 
Only about 15 per cent of these sources are likely to be 
real, equivalent to less than 0.5 per cent of the total. We 
therefore removed all possible side-lobes from our version of 
the catalogue prior to our analysis. 



2.2 The APM Catalogue 

The APM catalogue is derived from the UK Schmidt and 
POSS plates. For our present purposes we used APM data 
obtained from the POSS for the region covered by FIRST. 
The APM catalogue contains data from blue O plates and 
red E plates to limiting magnitudes of 21.5 and 20 respec- 
tively, with a positional accuracy of ~l". Each detection is 
parameterised in terms of position, magnitude, and a clas- 
sification into a stellar or a non-stellar object based on the 
light profiles of a detection; a 'merged' classification exists 
where an object falls between these two types. A more com- 
plete de scription of the APM sky catalogues ca n be found 
at URL ittp: //www. ast . cam. ac .uk/~apmcat/. 



2.3 Selection of Primary Sources 

In order to look for possible cases of multiple imaging we 
must define a parent sample of "primary" sources around 
which we will search for "secondary" images. 

With an average beam size of ~5", the FIRST sur- 
vey can confidently identify resolved structure down to an 
angular scale of ~2'.'5 (FWHM). For redshifts z ^ 0.2, 
this corresponds to a linear scale of J5 7 kpc (for Hq = 
65 kms~^Mpc~^). A region of this size will not vary in 
spectral or luminosity characteristics over periods of hun- 
dreds of years. We therefore placed a minimum angular size 
limit of 2" 5 along the deconvolved major axis of our primary 
sources. We chose the maximum angular size of our primary 
sources to be 15" to avoid confusion with potentially-lensed 
secondary images with similar angular separations. 

With the resolution of the FIRST survey many faint ra- 
dio sources will be resolved into multiple components with 
angular separations less than 15"; these components are 
listed separately in the FIRST catalogue. We therefore con- 
structed two classifications of primary sources, those with 
a single component and those with two components (see 
Figure hh. One could include primary sources with three or 
more components, but for simplicity we ignored such sources 
(inclusion of these sources would only increase our sample 
size by ~10 per cent). In the case of a single component our 
maximum primary source size of 15"refers to the major axis 
given in the FIRST catalogue, in the case of two components 
it refers to the angular separation of the components. We de- 
fined the central position of a two-component primary source 
to lie on the line connecting the individual component posi- 
tions, weighted by the integrated component fiux densities 
(see for example Macklin 1981). The "fractional centroid dis- 
tance" being the fractional distance to the centroid from the 
brighter component is defined as /3 = Fb/{Fa + Fb) where 



Fa , Fb are the integrated fiux densities of the brighter and 
dimmer components respectively. We can then write the an- 
gular distance from the brighter component to the weighted 
central position as, 6 x j3, where 6 is the angular separation 
of the two components. Two extreme cases of composite pri- 
mary sources are illustrated in Figure eI 

In order to select "clean" primary fields and to avoid 
double counting, we placed a window with a radius of 
15"centred on the central position (defined above) of each 
primary source. We only included a source in our primary 
sample if no additional components were within the annulus. 
This will help prevent more complex (i.e. 3-|- component) 
sources from infiltrating our sample. We finally imposed the 
condition that the integrated fiux density must be >35 mjy. 

To increase the probability of lensing, it is desirable to 
select background sources with as high a mean redshift as 
possible. To this end we selected only those primary radio 
sources with stellar optical counterparts in the APM cata- 
logue on the O and/or E plate. By choosing stellar objects 
with associated radio emission we expect to select a high pro- 
portion of q uasar s, and hence to maximise the mean redshift 
(see section 2.4.3). Radio sources with non-stellar i.e. galaxy 
identifications will tend to have redshifts ^0.2 and therefore 
have a greatly reduced probability of being lensed. Optical 
identifications with a limiting E magnitude brighter than 20 
also make it possible to characterise the redshift distribu- 
tion of the primary sample with optical spectroscopy using 
moderate-sized telescopes. In addition, an optical identifica- 
tion on a primary source will provide a filter against selecting 
clas sical double radio sources as lens candidates (see section 

A significant number of faint red galaxies are identified 
as non-stellar on the E plates, but have a stellar classifica- 
tion on the O plates due to the detection of only the central 
region of the galaxy. Figure illustrates the three distinct 
populations present in a primary source sample when we al- 
lowed there to be a non-stellar identification on one of the 
POSS plates. To simplify the optical identification process 
for ARCS, an optical object with different classifications on 
the two plates is not acceptable as a quasar identification. 
In fact, as Figure ra shows, those sources with stellar classi- 
fications on one plate and stellar or no classification on the 
other plate are predominantly blue, as would be expected 
for quasars (see for example the catalogue derived from four 
SDSS fihers described by Newberg et al. 1999). 

To take into account the positional accuracies of the 
APM catalogue, a 3" difference betw een th e APM and radio 
positions was tolerated (see section 2A3). Radio positions 
for the single component sources were taken directly from 
the FIRST catalogue, whereas the composite source posi- 
tions were as defined above. After applying all these selec- 
tion criteria, a final parent sample of 1023 primary sources 
remained (see Table 0) . 



2.4 Statistics of the Primary Source Sample 

2.4.1 Random Radio Associations Posing as 
Two-component Primaries 

By selecting two-component primary sources from a cata- 
logue with the depth of FIRST we will inevitably introduce 
a contamination to the primary sample in the form of pairs 



© 2000 RAS, MNRAS 000, 0-|l| 



4 P. M. Phillips et al. 



22 



20 



0) 
X3 



m li 



14 



g%. 



«o 



5 o 
< 

A 



; &- ++ 



Wo' ' ' 



^9 ^ 






+++ 




+A+ ! 



#^ 



f^t-^ V A+ 



16 



20 



17 18 19 

R magnitude 

Figure 3. Plot of R magnitude against B magnitude for a primary source sample where the APM optical classification must be stellar 
on at least one of the two plates. The diamonds (a) are sources with non-stellar classification on the E (red) plate, the triangles (b) are 
sources with non-stellar classification on the O (blue) plate and the crosses (c) are sources with a stellar classification on both plates. The 
range of the plot is restricted to illustrate the main concentration of sources (~65.5 per cent of all those with optical identification) - most 
of those not plotted have a magnitude measured on one plate only. Three separate populations are clearly visible: (a) (~15.5 per cent) 
contains relatively red sources which are probably faint red galaxies; (b) (~6 per cent) contains a small population of generally blue 
sources; (c) (~78.5 per cent) are relatively blue stellar sources with a distribution that would be expected for a population of QSOs (for 
example Newberg et al. 1999). For ARCS, only objects in category (c) were selected (see text). 




Table 1. Synopsis of the source selection stages of the ARCS 
project. 



2 3 

Separation (arcsec) 

Figure 4. Histogram of the separation between the optical and 
the radio position for all primary sources. The dashed line shows 
the increasing number of random associations with search radius, 
while the solid line shows the best fit Rayleigh-plus-random dis- 
tribution, for which a-=0'.'86. Whilst our selection criterion allows 
a maximum optical-radio separation of 3", we have extended this 
plot to 5" to illustrate the increasing dominance of random asso- 
ciations above 3". 



FIRST catalogue, 98feb04 

After radio selection (1/2-component) 

After optical selection (1/2-component) 

Total ARCS primary sample 

Initial lens candidates 

VLA follow-up candidates (after 

by-eye examination of FIRST maps) 

MERLIN follow-up candidates 

Number of lens systems 



Number of sources 

382 892 

7668/5264 

627/396 

1023 

91 

38 
2 




of sources which have no physical association. We expect this 
number to be small, but we have made a rough estimate of 
the number of chance associations of unrelated radio sources 
posing as two-component sources in the FIRST catalogue. 
Our selection routine found 5264 two-component primary 
sources from the radio selection criteria alone i.e. prior to 
the application of the optical criteria. From the known sur- 
face density of FIRST sources the number of random asso- 
ciations was estimated to be ~100. This is a contamination 
of only ~2 per cent amongst the two-component sources be- 
fore our optical criteria have been applied (i.e. <1 per cent 
of all primaries - see Table y). Within our final primary 
sample we expect this percentage to have been reduced sig- 
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nificantly because randomly associated pairs of sources will 
be less likely to have stellar counterparts at their centroid 
position in comparison to a genuine physical pair. Our con- 
clusion is that the level of contamination to the primary 
sample from the selection of false two-component primary 
sources will be negligible. 



2.4-2 Random Radio - Optical Associations 

Since the error in the radio astrometry is negligible relative 
to the optical astrometric accuracy, and we assume that the 
optical positional error is purely gaussian, then the separa- 
tion between the optical position and the true (i.e. radio) 
position, r, is described by a Rayleigh distribution, viz. 



P(r)dr = ^exp 
a'' 



r 



dr 



where a is both the mean and the standard deviation of 
the distribution. However, there is a small but significant 
contamination within our primary sample of radio sources 
with chance associations with, for example, galactic stellar 
sources. This modifies the number of optical 'identifications' 
by 

Nran{r)dr = Ei.Nrad'iTVrdr 

where S* is the stellar surface density arcsec"'^ in the 
APM catalogue from both POSS plates, N^ad is the num- 
ber of primary sources before optical identification, and r 
is the distance from the radio source. Given the fact that 
Et ~ 5 X 10""' arcsec"^ within the FIRST region, the total 
number of random associations amongst the 1023 primary 
sources is estimated to be ~180 to a radius of 3". Figure 
shows a histogram of the angular separation of optical 
and radio positions in our primary source sample. The dis- 
tribution is well described by a Rayleigh-plus-random dis- 
tribution, increasing our confidence that we are selecting 
extragalactic radio sources with quasar identifications. The 
best fit value of a =0'.'86 agrees well with the quoted error 
in the APM positions. However, to select as many primary 
sources as possible we retained our 3" tolerance for FIRST- 
APM positional discrepancies. In an analysis of the lensing 
frequency within the sample, we must modify the primary 
sample size to ~850 sources to take account of the expected 
number of random associations. 



2.4.3 Mean Redshift 

We searched the NASA/IPAC Extragalactic Database 
(NED) for additional information on our primary sources. A 
total of 133 sources were found to have a measured redshift 
at the FIRST position. We calculated values of z ~ 1.09, 
R ~ 17.5 and B ~ 18.1 for this sub-sample. The total 
primary sample is on average fainter with R ~ 18.6 and 
B ~ 19.5. Whilst the sub-sample obtained from NED will 
be unquantifiably contaminated by a multitude of selection 
effects, it is reasonable to assume that the total primary 
sample will have a higher value for z. In fact, an initial anal- 
ysis of spectroscopic data from a representative but limited 
sub-sample of the primary sample suggests that z ~ 1.4 
(Phillips et al., in preparation). 



Primary source 



Secondary source 
> 20 per cent integrated 
flux density of primary 




2".5-1.5" 




Figure 5. Illustration of an ARCS lens candidate. Note that 
while the secondary source in this case is depicted as being smaller 
than the primary source, in practice no size criteria were applied 
to the secondary sources. 



2.5 Initial Lens Candidate Selection 

Having selected the primary sources the resolution and sensi- 
tivity of the FIRST survey allowed us to search for potential 
multiple imaging events directly. We looked in the catalogue 
near each primary source for any secondary sources within 
an angular separation range of 15" and 60" (Figure Bl); any 
such secondary sources are putative lensed images of the 
primary sources. We did not impose any size selection cri- 
teria to the secondaries because a secondary lensed image, 
with a lower magnification than the primary, will not nec- 
essarily be resolved by the FIRST beam. We also omitted 
optical selection criteria; since the secondary will often be 
weaker, its optical emission could place it below the limit 
of the POSS/APM catalogue; in addition a stellar object 
may have varied in magnitude. However, we did impose a 
maximum fiux ratio of 5:1 between the primary and the sec- 
ondary. This places a lower limit on the integrated fiux den- 
sity of a secondary source in the weakest systems to 7 mjy. 
If a secondary source is indeed a lensed image of the pri- 
mary, the conservation of surface brightness means that the 
secondary will not be resolved out i.e. it will be detectable 
in FIRST. 



In section 2.3 we mentioned that a large proportion 



of extragalactic radio sources will be resolved into multiple 
components by the FIRST survey. Many such sources consist 
of two 'lobes' of emission, the so-called "classical doubles" . 
It is possible that in selecting primary sources with nearby 
secondaries, we are selecting a high proportion of classical 
double radio sources as lens candidates. However, the opti- 
cal identification criterion applied during the primary source 
selection process provides a strong filter against these types 
of radio source posing as lens candidates. This is because 
the optical counterpart to such an object would be expected 
to lie between the lobes rather than on one or both of them 
as our selection criterion demands. 

These selection criteria produced only 91 initial lens 
candidates i.e. over 90 per cent of the primary sample could 
be immediately rejected using the FIRST and POSS/APM 
information. 
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Figure 6. Examples of FIRST maps from our 91 initial lens candidates. Stellar objects are shown as crosses with a ±3" extent, to 
illustrate the optical identification tolerance; non-stellar objects are shown as diamonds whose dimensions describe the fitted ellipse 
parameters in the APM catalogue. Radio contours are logarithmic with base 2 starting at Sit. Labels P and S denote primary and 
secondary components respectively. The sources are 0713+369 (left) and 1358+290 (right). The primary 'source' is located at the centre 
of the map in each instance. It is clear that 0713+369 and 1358+290 are single radio sources: 0713+369 has bridge emission between 
the primary (central) component and the two secondary components are edge brightened away from each other; similarly 1358+290 has 
a central primary component with a southern secondary component edge brightened away from the primary. The northern secondary 
component is likely to be a counter-lobe. Further examples of FIRST maps of ARCS sources are shown in Figures m, H and U 



2.6 Further Rejection of Lens Candidates 

Up until this point the lens candidate selection process was 
conducted automatically. We then went on to examine the 
FIRST and APM data individually by eye. 

The FIRST radio images of the 91 initial lens candidates 
along with APM information for each field were combined 
using the AIPS (Astronomical Image Processing Software) 
package. Examples of typical fields are shown in Figure H and 
the top of Figures [?|, g and |9|. It is immediately clear that 
we can reject some more candidates by eye. The arguments 
for rejecting a system as a lens candidate can be broken 
down into two categories: (1) those where the primary and 
secondary sources are clearly components of a single larger 
source; and (2) those where lensing arguments can be in- 
voked. 

(1) Single source arguments: 

The single source rejection arguments rely on the nature of 
known multiple-component radio sources. The morphologies 
of such sources can include specific features which would not 
be expected if the components were lensed images of each 
other. Examples are: extended emission bridging the region 
between the primary and secondary source, and 'outer' edges 
of a primary-secondary pair being brighter than the 'inner' 
edges. We also reject lens candidates where we find an opti- 
cal counterpart between the primary-secondary pair coinci- 
dent with an additional unresolved radio component; this is 
the position where we would expect an AGN to be present. 
Examples of sources rejected as single sources on the basis 
of the FIRST maps are 0713+369 and 1358+290 shown in 
Figure |6[ 



(2) Lens arguments: 

We use the basic fact that the surface brightness of the back- 
ground source is conserved in each lensed image. However, 
we can only apply this argument to the extended radio struc- 
ture. The image fiux ratio in the radio and optical regimes 
can be differe nt d ue to the the variability arguments laid 



out in section 2.5 



From the original 91 sources with secondary compo- 
nents, 53 could be rejected following this individual exami- 
nation leaving only 38 lens candidates which required further 
observations i.e. less than 4 per cent of the primary sample. 



3 RADIO FOLLOW-UP OF CANDIDATES 

In order to test the lens hypothesis for the remaining candi- 
dates we made higher resolution observations to study can- 
didate source morphology and surface brightness in more 
detail. 



3.1 Stage 1 — VLA Observations 

The sample of 38 lens candidates was observed with the VLA 
in A configuration at 6 cm and 20 cm in two 12-hour periods 
on 1 and 7 April 1998. The observing bandwidth was 50 MHz 
in each of two IF channels for both wavelengths, with ob- 
serving periods of approximately 8 minutes per source per 
wavelength. The rms sensitivity was ~ 0.05 mjy at both 
wavelengths. The minimum beam sizes obtained were 0'.'4 
at 6 cm and l'.'4 at 20 cm when sources were observed near 
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Figure 7. Examples of single (i.e. intrinsic double) sources from the ARCS VLA follow-up sample: 1248+393 (left) and 1519+345 (right). 
From top to bottom: FIRST radio map; VLA A configuration 20 cm; VLA A configuration 6 cm. Contours and symbols are as Figure pi 
In both cases it is unlikely that our optical identifications are correct because they are coincident with 'lobe' emission. A candidate for 
the true identification of 1248+393 is located between the primary and secondary radio sources. 
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Figure 8. Examples of pairs of unrelated sources from the ARCS VLA follow-up sample: 1000+417 (left) and 1034+379 (right). From 
top to bottom: FIRST radio map; VLA A configuration 20 cm; VLA A configuration 6 cm. Contours and symbols are as Figure pi 
In both examples the position of the optical identification of the primary source is consistent with the radio morphology. In addition 
1034+379 exhibits a non-stellar optical counterpart between the two lobes of the secondary source. 
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Figure 9. The two remaining lens candidates from tlie ARCS VLA follow-up sample: 1053+457 (left) and 1519+387 (right). From top 
to bottom: FIRST radio map; VLA A configuration 20 cm; VLA A configuration 6 cm. Contours and symbols are as Figure H. In both 
examples the position of the optical identification of the primary source is consistent with the radio morphology. 
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Figure 10. High resolution 6 cm MERLIN maps of the final lens candidates 1053+457 (top) and 1519+387 (bottom). Individual maps of 
the secondary (left) and primary sources (right) are shown. In the case of the composite primary source of 1519+387, only the unresolved 
component is shown. Contours are logarithmic with base \/2 starting at 3a-. 



the zenith. The data were calibrated in a standard fash- 
ion using the AIPS package. Flux calibration was based on 
two 2-minute observations of 3C286 assuming a flux den- 
sity of 7.360 Jy at 6 cm, and 14.554 Jy at 20 cm based 
on the expression derived by Baars et al. (1977). Polarisa- 
tion calibration was based on three 1.5-minute observations 
of 1613+342, which was assumed to be unpolarised to cor- 
rect for antenna effects, while 3C286 was used for phase 
offset corrections between the two measured orthogonal po- 
larisation systems. Each source was CLEANed at both 6 cm 
and 20 cm for Stokes parameter I using the DIFMAP soft- 
ware package (Shepherd 1997) and where appropriate self- 
calibration was applied. Polarisation and spectral index data 
reduction was performed with the standard AIPS routines. 
In these cases a uniform convolving beam size of l'.'5 was 
used at both 6 cm and 20 cm for ease of image comparison. 
A Gaussian taper was applied in the uv plane to the 6 cm 



data to 30 per cent at 200 kA in order to restrict the uv 
range to that of the 20 cm data. 

We used similar arguments to those described in sec- 
tion 2.6 to classify each system. The 38 lens candidates 



fall into three distinct groups: single 'classical double' radio 
sources (~68.5 per cent); unrelated small angular separa- 
tion pairs (~26.5 per cent); and 2 remaining lens candidates 
(~5 per cent). Figures M,m and H show examples from each 
of these groups. 

Examples of single sources (Figure i7|): 

In both examples, 1248+393 and 1519+345, it can be seen 
from the A configuration 20 cm data that the primary and 
secondary components are outer edge-brightened. The 6 cm 
data show that the spectra of the outer edges are flatten- 
ing, consistent with these regions being hot-spots embedded 
within diffuse lobes. 
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Examples of two unrelated double sources (Figure H): 
The 6 cm A configuration map of example 1000+417 shows 
a compact double primary (northerly) source with a cen- 
tral optical counterpart. The secondary source consists of 
a central unresolved flat-spectrum component and extended 
components to the north-east and south-west, with extended 
emission bridging the entire source. We can rule this example 
out as a lens candidate using surface brightness arguments. 
In example 1034-1-379 the 20 cm A configuration data 
shows extended emission bridging two components in the 
primary (northerly) source, while the 6 cm data reveals an 
unresolved flat-spectrum component at the position of the 
optical counterpart. The secondary source at 20 cm contains 
two edge-brightened components with a non-stellar optical 
counterpart between, in a position that a host galaxy would 
be expected. This example can be ruled out as a lens can- 
didate based on the dissimilar morphologies of the primary 
and secondary sources, and the orientations of each source 
relative to one another. 

Remaining lens candidates (Figure H): 

The A configuration map of 1053-1-457 at 20 cm shows two 
compact double component sources in a plausible lens con- 
figuration. The 6 cm data show that the primary (west- 
ern) source in fact consists of extended emission bridging 
two unresolved components, whereas the secondary (east- 
ern) source has a single unresolved component with an ex- 
tended component to the south-east. Despite the dissimi- 
larities between sources, we cannot unequivocally rule out 
this system as a lens candidate since unresolved radio com- 
ponents can vary significantly over the time-scale of a time 
delay. The 1519-1-387 primary (western) source is a com- 
posite consisting of an extended western component and an 
unresolved eastern component in the A configuration 20 cm 
data. The secondary source appears largely unresolved in the 
6 cm data with a possible detection of extended emission 
towards the unresolved component of the primary source. 
While these components could be interpreted as a single 
source with a central unresolved core at PI (see Figure bl), 
high resolution MERLIN data are needed to confirm the ex- 
tension in the secondary source. At this stage it remains as 
a lens candidate. 



3.2 Stage 2 - MERLIN Observations 

High resolution data on the two remaining lens candidates 
1053+457 and 1519+387 were taken with MERLIN at 6 cm 
in two 16- hour long-track observations on 12 and 21 De- 
cember 1998. The data were in 16x1 MHz bands with an 
integration time of 8 seconds, giving an rms sensitivity of 
~ 0.05 mjy; the minimum beam size was ~0'.'04. Flux cal- 
ibration was based on 25 minute observations of 3C286 
assuming a flux density of 7.09 Jy. Initial editing of the 
data was carried out with the standard MERLIN programs 
with calibration and mapping again performed in AIPS and 
DIFMAR 

The unresolved components in the VLA 6 cm data 
of the primary source of 1053+457 are both clearly re- 
solved with MERLIN (Figure hol) and appear to be edge- 
brightened. Additionally a central unresolved component is 
revealed between the two extended components and is close 



to an optical object (see Figure H). The unresolved com- 
ponent in the secondary source also appears to be slightly 
resolved while the extended component to the south east is 
almost resolved out. On the lens hypothesis we would expect 
the two resolved components in the primary source also to 
be detected in the secondary source. This is not the case. 
We detect only one resolved component in the secondary, 
which is also of lower surface brightness than either of the 
resolved components in the primary source. We can therefore 
reject 1053+457 as a lens candidate. The secondary source 
in 1519+387 is clearly resolved with MERLIN while the pri- 
mary component remains unresolved. We can categorically 
reject this source as a lens candidate. 



4 DISCUSSION AND SUMMARY 

We have completed what is currently the largest unbiased 
survey for multiple image gravitational lens systems with im- 
age separations >15". Our primary sample of 1023 sources 
was selected from the FIRST and APM catalogues, and our 
search for lens candidates within this sample was achieved 
by looking for additional nearby images within the FIRST 
catalogue. As a result our follow-up required only 38 can- 
didates be observed with the VLA, and just two of these 
necessitated higher resolution observations with MERLIN. 
No cases of gravitational multiple imaging with image sepa- 
rations in the range 15"to 60"were found. A synopsis of the 
numbers of sources at each stage of the selection process is 
given in Table nl 

Wambsganss et al (1995) made predictions of "wide- 
separation" lensing rates using A'^-body simulations with 
SCDM parameters (see also Wambsganss, Cen & Ostriker, 
1998). However, the resolution and scale of the simulations 
were not adequate to provide a reliable estimate of the num- 
ber of multiple imaging events that we would expect to see in 
an unbiased search. In particular, cluster mass density pro- 
files must be simulated more accurately because the shape 
of the profile has a significant influence on a cluster's ability 
to produce multiple images. In order to draw any significant 
cosmological conclusions from our result more realistic pre- 
dictions of the number of wide-separation multiple imaging 
events in different cosmological scenarios need to be made. 
The results of such simulations are also necessary if we are 
to determine the sample size for any future expansion of the 
ARCS project. 

To examine the practicability of increasing the size of 
the ARCS sample we have calculated the maximum size 
of a sample which can be selected from currently available 
FIRST data and the pending SDSS data. As described in 
section 2.1 FIRST has been designed to coincide with the 
SDSS survey region. It is estimated that ~50 per cent of all 
sources in FIRST wiU be identifiable on SDSS to m{v) ~24 



down to the 1 mJy limit of FIRST. Using our radio-only 
primary source sample of 12 932 (Table [D) this could result 
in a primary sample of >6500 if SDSS is used instead of the 
APM catalogue at the optical identification stage. Taking a 
more conservative approach by specifying an optical mag- 
nitude limit of miv) ~23 we would still expect to increase 
our primary sample to >3000. It should be noted that these 
figures are firm lower limits on what could be achieved be- 
cause our radio-only primary sample has a lower 1.4 GHz 
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flux limit of 35 mjy (c.f. the 1 mjy limit of FIRST). The 
latest version of the FIRST catalogue dated 99jul21 con- 
tains 477 247 reliable sources over 5450 square degrees of the 
northern sky (c.f. 382 892 sources in the 98feb04 version we 
used) . If we also lower the integrated flux density limit of a 
primary source to 5 mJy, thus giving a minimum secondary 
flux density at the detection limit of the FIRST survey i.e. 
1 mJy, a total radio-only primary sample of 62 891 sources 
results. Applying the above optical arguments with a limit- 
ing magnitude of m(v) ~23, this should produce a primary 
sample of > 15 000. It would appear that upon completion 
of SDSS there will be no problem constructing a sample of 
background sources at least an order of magnitude larger 
than that presented in this paper. 

It is practicable to search such a primary source sam- 
ple for instances of wide-separation lensing events. With the 
experience gained in this pilot study, we believe that the ef- 
flciency of our search technique could be improved by up to 
a factor ~2. This is largely based upon our increased confi- 
dence in the FIRST survey images for use in rejecting false 
positives before any further observations of lens candidates 
are made at higher resolution. Selecting a primary sample 
of, say, 10 000 sources, a two-fold improvement in the initial 
candidate selection would result in only ~200 candidates re- 
quiring VLA A configuration follow-up with ^--^10 requiring 
further observations with MERLIN. This is equivalent to ~3 
days of VLA time and ~3 days of MERLIN time. 

In summary, we have established that it is possible to 
make a reliable unbiased search for lens systems with image 
separations in the range 15"to 60". Furthermore it is pos- 
sible to use existing radio and optical data to significantly 
reduce the required observing time for such a search. What 
we currently lack are reliable predictions of the lensing fre- 
quency for different cosmologies to compare with our null 
result. 
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